Coexistence of the antiferromagnetic and superconducting order and its effect on spin 

dynamics in electron-doped high-T c cuprates 
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In the framework of the slave-boson approach to the t — t' — t" — J model, it is found that for 
electron-doped high-T c cuprates, the staggered antiferromagnetic (AF) order coexists with super- 
conducting (SC) order in a wide doping level ranged from underdoped to nearly optimal doping at 
the mean-field level. In the coexisting phase, it is revealed that the spin response is commensurate 
in a substantial frequency range below a crossover frequency ui c for all dopings considered, and it 
switches to the incommensurate structure when the frequency is higher than lo c . This result is 
in agreement with the experimental measurements. Comparison of the spin response between the 
coexisting phase and the pure SC phase with a d x 2_ y 2-w&ve pairing plus a higher harmonics term 
(DP+HH) suggests that the inclusion of the two-band effect is important to consistently account for 
both the dispersion of the spin response and the non-monotonic gap behavior in the electron-doped 
cuprates. 

PACS numbers: 74.20.Mn, 74.25.Ha, 75.40.Gb 



I. INTRODUCTION 



The pairing symmetry of the hole-doped high-T c su- 
perconductors is generally believed to have the dominant 
rf a .2_j / 2-wave pairing. However, the pairing symmetry of 
the electron-doped high-T c superconductors is still under 
debateii»2iMiSi2i2iS While no consensus has been reached 
yet, more and more recent experimental results have 
suggested that the order parameter of electron-doped 
cuprates is likely to have a dominant d x 2_ y 2-w&ve pairing 
symmetry— diSiSiZiS and with an unusual non-monotonic 
gap function. 

Although various explanations have been pro- 
posed to account for the non-monotonic behavior, 
they can generally be categorized into two scenar- 

^,5 ^ 0,11 ^ 13,14,15,16,17,18 Qne ^ tQ extend ^ gu _ 

perconducting (SC) gap out of the simple d 2 .2__„2-wave 
via the inclusion of higher harmonics termi^i 7 -^ ' 1 ^ 11 ' 12 ! 13 
From theoretical perspective, the non-monotonic d x 2_ y 2- 
wave gap appears under the assumption that the d x 2_ y 2- 
wave pairing is caused by the interaction with the contin- 
uum of overdamped antiferromagnetic (AF) spin fluctu- 
ations. In this scenario, the non-monotonic gap behavior 
is described by the combination effect of a d 2 .2_ y 2-wave 
paring plus a higher harmonics term (DP+HH). There- 
fore, it is an intrinsic property of the SC state regardless 
of the presence of the AF order, and a simple one-band 
model is used to reproduce the non-monotonic gap be- 
havior. The other argues that the non-monotonic be- 
havior is the outcome of the coexistence of the AF and 
the SC ordersi 14 i 15 ' 16 i 17 ' 18 This scenario assumes that the 
AF order disguises the d x 2_ y 2-w&ve character of SC gap. 
When the AF order is introduced, the resulting quasi- 
particle (QP) excitation can be gapped by both orders 
and behaves to be non-monotonic, although the SC gap 
itself is monotonic. The scenario gained support from 
the angle-resolved photoemission spectra (ARPES) mea- 



surements, where two inequivalent Fermi pockets around 
(tt,0) and (ir/2,n/2) have been detectedJ^ This phe- 
nomena is well explained in terms of the k-dependent 
band-folding effect associated with an AF order which 
splits the band into upper and lower branches , 14 i 20 i 21 i 22 
leading to the two-band and/or two-gap model. 

Recently, neutron scattering experiments in electron- 
doped cuprates have revealed that the spin response is 
commensurate in a substantial frequency range below 
a crossover frequency u> c ^i 24 ! 25 , 26 ! 27 , 28 which consti- 
tutes a distinct difference from the widely studied hour- 
glass dispersion in the hole-doped cuprates^ Although, 
both scenarios mentioned above can account for the non- 
monotonic gap behavior of the electron-doped cuprates, a 
comparative study on the spin dynamics between the two 
scenarios is deserved to demonstrate the possible differ- 
ences and therefore serve to select the reasonable model 
for electron-doped high-T c cuprates. 

In this paper, we investigate the spin dynamics in the 
coexisting phase of the AF and the d x 2_ y 2-w&ve SC or- 
ders. The calculation is based on a self-consistent deter- 
mination of the QP dispersion, the AF order and the SC 
gap at the slave-boson mean-field level of the t—t'—t" — J 
model. It is shown that the AF and SC orders compete 
and coexist in a substantial doping range in the under- 
doped regime. The spin response is commensurate below 
a crossover frequency lu c for all dopings considered, and 
it becomes incommensurate when the frequency is higher 
than oj c . This result is qualitatively consistent with ex- 
perimentsj 23 i 24 ' 25 i 26 ' 27 i 28 While in the framework of the 
pure SC state with d x 2_ y 2-w&ve and/or DP+HH, ^LSLSL, 
though an extended region of a commensurate spin fluc- 
tuation also exists, it evolves into an incommensurate 
spin fluctuation at low frequencies, which is not consis- 
tent experiments. Therefore, our result suggests that the 
inclusion of the two-band effect resulting from the coex- 
isting AF and SC orders is important to consistently ac- 
count for both the spin dynamics and the non-monotonic 
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gap behavior in the electron-doped cuprates. 

The paper is organized as follows. In Sec. II, we intro- 
duce the theoretical model and carry out the analytical 
calculations. In Sec. Ill, we present the numerical results 
with some discussions. Finally, we present the conclusion 
in Sec. IV. 



J 



II. THEORETICAL MODEL 



The Hamiltonian of the two dimensional t — t' — t"- 
model on a square lattice is written in the form, 



J 



H = -t J2 {c\ a c ja + H.c.)-t l (clc ja + H.c.)-t 2 KW + ff-c) 

<ij>,a <ij>2-A <ij>3,& 

+J ( S * ' S i - J n i n i) - Mo c\ a d<y (1) 



<i]> 



Where the summations < ij >, < ij >2, < ij >3 
run over the nearest-neighbor(n-n), the next-n-n, and the 
third-n-n pairs respectively, Si is the spin on site i. This 
Hamiltonian can be used to model both hole-doped and 
electron-doped systems after a particle-hole transforma- 
tion. For electron doping, one has t<0, ti>0 and t2<0. 
The slave-boson mean-field theory (SBMFT) is used to 
decouple the electron operators Ci a to bosons bi carry- 
ing the charge and fermions fi a representing the spin. 

Then, the local constraint b\bi + J2i& ftrfi<r = ^ * s satis- 
fied averagely at the mean-field (MF) level. We choose 
the spinon pairing order Ay=< / iT /^ - / a / iT >=±A 
, where Ay = A(— A) for bond < ij > along x{y) di- 
rection, the uniform bond order Xij = J2a < furfjo >= Xi 
the AF order < f^fy - fj^ > /2 = (-l) l rn, and re- 
place bi by < bi >=^fx due to boson condensation. After 
the Fourier transformation, the mean-field (MF) Hamil- 
tonian can be written in the Nambu representation, 

H = H ° 1 OtiMtyCQ*) + 2N J(X 2 + m 2 + A 2 /2) - Nfi, 

k 

(2) 

here, the Nambu operator C'(k) = 
(fk-\> fl+Qp f-*l> f-*-Ql)> and 
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-2Jm 


-JA k 
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JA k 


JA k 
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-2Jm 





JA k 


-2Jm 


- £ k+Q 



(3) 



where, e k = (— 2tx — Jx)(cosk x + cosk y ) — 
At\x cos k x cos k y — 2tix{co% 2k x + cos 2k y ) — /i and A k = 
A (cos k x — cos k v ). fj, is the renormalized chemical poten- 
tial, N is the total number of lattice sites, and Q = (ir, ir) 
is the AF momentum. Note that the wave vector k is 
restricted to the magnetic Brillouin zone (MBZ) in all 
follows. 

Diagonalizing of the Hamiltonian (2) by an unitary 
matrix U(k) leads to four energy bands -E'i(k) = E£ , 



E 2 (k) = E k , E 3 (k) = -E k , E 4 {k) = -E+, with 

E£ = y/(g) 2 + (JA k )*, (4) 

where g = e k ± \/(e k ) 2 + 4J 2 m 2 with e± = (e k ± 
£k+cj)/2. And the free energy is written down (Boltz- 
mann constant ks = 1), 

E v 

F = -2T Y ln[2cosh(^)]-^iV+2iVJ(x 2 +m 2 +A 2 /2). 

k,v=± 

(5) 

The MF order parameters %, A, m and the chemical po- 
tential /i for different dopings x can be calculated from 
the self-consistent equations obtained by dF/d\ = 0, 
dF/dA = 0, dF/dm = 0, and dF/dfi = -N(l - x), 
respectively. The magnitudes of the parameters are 
chosen as t=-3.QJ, ti=1.02J, t 2 =-0.51J and J=I00 
meV to model the Fermi surface observed in ARPES 
experiment j 19 i 20 

Then, the bare spin susceptibility (transverse) is given 
by, 

xo (q>q>) = ^ < S+(r)Sl q ,(0) > (0) , (6) 

where < • ■ • >( ) means thermal average over the eigen- 

states of H, 5+ = Ylk fk+qff^l ^ s * ne spin operator. 
Considering that k is restricted to the MBZ, an explicit 
calculation shows that the spin susceptibility should be 
expressed in the following matrix form, 

i±( ne ,\-{ x^(q,q,^) xJ(q,q + Q,^) \ (7] 
Xolq ' >~ Uo(q + Q.q,^) x ± (q + Q,q + Q^)y A ' 

where, the nondiagonal correlation function xt with q' = 
q + Q arises due to the umklapp process. The matrix 
elements of the bare spin susceptibility, which come from 
the particle-hole (p — h) excitations, are given by, 
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. 2 2 



/(£ m (k))-/(S„(k + q)) 



ij — l m,n— 1 k 

i-/(E m (k))-/(£;, 



£„(k + q) 
k + q)) 



+ a 2 - 



/(£ n (k + q))-/(£ m (k)) 

i w- J B„(k + q ) + i; m (k) + ir 



w - 



£;„(k + q ) + s m (k) + ir 



E rn {k)+iT w-£„(k + q)+£ m (k) 
, /(g m (k)) + /(£? n (k + q))-l 

2 w-£„(k + q )-i; m (k) + ir J 



(8) 



where, /(-Ek) is the Fermi function and 



oi = ^n( k + q)^0+r,'-,,)«(k + q) U lm (k) J7* m (k) + [/* n (k + q)t/ (i+2 )„(k + q)f/ im (k)[/ ( * +2+r) , _ v)m (k) , 

«2 = f/(i+2)n(k + q)^(*,+2+r,'-,,)n( k + q) +2)m ( fc ) C/ (i+2)m ( k ) + ^(i+2)n( k + q)C/* n (k + q) C/ ( * +2)m (k) L/ (j+r; , _„) m (k) , 

fo l = ^m( k + q)^0+ J ?'- J ?)n( k + q)t / (*+2)m( k ) ?7 0'+2)m(k) - U* n (k + q) U(j +2)n (k + q) C/ ( * +2)m (k) [7 0+t? , _„ )m (k) , 

62 = C/ (4 +2)«( k + q)U{ ]+2+ , n ,_, n)n (k + q)U im (k)U* m (k) - U (l+2)n (k + q)C/;„(k + q)C/ jm (k)[/ ( * +2+t? ,_ t?)m (k). (9) 



The renormalized spin susceptibility due to the spin 
fluctuations is obtained via the random-phase approxi- 
mation (RPA), 

X ± (q,-) = T -%^- ! (10) 

where, 

J - ( J(q) ° ^| (11) 
Jq -{ J(q + Q)J 

with J(q) = ./(cosg.,, + cosq y ). In the coexisting phase 
of the AF order and SC order, a is taken as 1. As for 
the pure SC state with DP+HH, we choose a slightly 
small a = 0.72, the criteria for choosing a is to set the 
AF instability at x = 0.12. The parameter T=0.04J is 
introduced to account for the QP damping rate which 
comes from the scattering off other fluctuations that are 
not included here. 

III. NUMERICAL RESULTS AND DISCUSSION 

In Fig. 1, we show the MF parameters m an d A as 
a function of doping x. For a comparison, we also show 
the doping dependence of the MF SC gap Ai obtained 
without considering the AF order by setting m — 0. It is 
seen that the staggered magnetization m decreases with 
increasing doping x, and goes sharply to zero at x « 0.16, 
which implies a phase transition from the antiferromag- 
netism (AFM) phase to the paramagnetic phase. The SC 
order parameter, on the other hand, increases its value 
initially up to an optimal doping level, and then decreases 
upon further doping, forming a generic SC dome^ 3 How- 
ever, the SC order parameter Ai without the inclusion of 



the AF order exhibits a monotonic decrease with doping, 
which deviates obviously from the experimental obser- 
vations. Furthermore, the SC order parameter A with 
an AF order shows a noticeable suppression compared to 
Ai, exhibits a competitive character with the AF order. 
But, they also coexist in a substantial doping range. The 
MF phase diagram also shows that the optimal doping 
is rather low compared to that deduced from the experi- 
ments. This may be due to the fact the SBMFT includes 
only the MF value of the order parameters and treats the 
no-double occupancy on the average. However, the sim- 
ilarity of the phase diagram obtained by the SBMFT to 
that of the variational quantum-cluster theor y 16 i 22 vali- 
dates the SBMFT as a low energy effective theory. Here 
our aim is to use the MF theory as an effective model 
to study the effect of the AF order on the spin dynam- 
ics, and then to compare the two-band and/or two-gap 
model with the simple one-band model. Therefore, the 
relatively simple SBMFT is qualitatively enough for our 
purpose. We note that a similar phase diagram has been 
obtained before*^ 

The doping dependence of the renormalized spin sus- 
ceptibility Imx(q, u>) at a low frequency u> — 0.04 J in 
the coexisting phase is presented in Fig. 2. In this fig- 
ure, it is found that the low-energy excitations exhibit 
commensurate peaks for all x, which consists with the ex- 
periments well 2 ^. The inset shows the spin susceptibility 
Im%(q, a;) at doping x = 0.15 in the pure SC state with 
DP+HH which is used to reproduce a non-monotonic SC 
gap behavior. One can see that the spin response is in- 
commensurate at low frequency without considering the 
AF order. 

Detailed frequency dependence of the spin response 
in the coexisting phase and the pure SC phase with 
DP+HH at doping x — 0.15 are shown in Figs. 3(a) 
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and 3(b), and Figs. 3(c) and 3(d), respectively. The dif- 
ference in the low frequency regime of the two phases is 
more evident here. The spin fluctuation is commensu- 
rate in a substantial frequency range below a crossover 
frequency tu c rs 0.52 J and down to the lowest frequency 
considered in the coexisting phase, and it switches to 
be incommensurate when the frequency is higher than 
u>c ss 0.52 J [Figs. 3(a) and 3(b)]. This feature agrees 
with the neutron-scattering measurements on electron- 
doped cuprates that have been reported recently^ While 
for the pure SC phase with DP+HH, the spin response 
is incommensurate at low frequency, then it switches 
to be commensurate within the intermediate frequency 
range, and becomes incommensurate again at higher fre- 
quency— These results can be summarized in the in- 
tensity plot of the imaginary part of the renormalized 
spin susceptibility Imx(q, w) as a function of frequency 
and momentum along (jr, q y ) direction, as shown in Fig. 

4. In the figure, the solid line indicating the peak posi- 
tion is the dispersion of spin excitations. The commen- 
surate spin fluctuation prevails below u> c for the coex- 
isting system [Fig. 4(a)]. For the pure SC phase with 
DP+HH, the dispersion shows a hourglass-like behavior 
[Fig. 4(c)], which is similar to the hole-doped one, and 
does not consistent with the experiments on electron- 
doped cuprates^ 

In the presence of the AF order, the energy band of 
quasiparticles is split into two bands. Therefore, the 
particle-hole excitations that contributed to the spin 
susceptibility are composed of two kinds of excitations, 
the intra-band and the inter-band excitations. In Fig. 

5, we present the results for the bare spin susceptibil- 
ity xo(q, ^) (without the RPA correction) coming from 
the intra-band and the inter-band contributions, respec- 
tively. Figs. 5(al) and 5(a2) denote the imaginary 
part of Xo(q>w), Figs. 5(bl) and 5(b2) the real part. 
One obvious feature is that, the intra-band contribu- 
tion is zero at the AF momentum Q = (tt, it), leading 
to the incommensurate spin response. It results from 
the fact that the coherence factor in the spin suscepti- 
bility due to the intra-band excitations, 1 — [(2Jm) 2 — 

Ek+qEkl/^k+q + (2Jm) 2 A/4 + (2Jm) 2 ] [where, e k = 
(—2tx — Jx)(cosk x + cos k y )] is zero at Q. While, the 
inter-band contribution is commensurate for all frequen- 
cies. At low frequencies, the inter-band excitations have 
a larger contribution to the spin susceptibility than the 
intra-band excitations, so the spin fluctuation is com- 
mensurate. However, with the increase of frequency, the 
intensity of Imxo (q, uj) due to the intra-band contribu- 
tions increases more rapidly than the inter-band contri- 
bution. As a result, the spin fluctuation switches from a 
commensurate to an incommensurate structure. 



IV. CONCLUSION 

In this paper, we have investigated the spin dynamics 
in the electron-doped cuprates in the coexisting phase of 
the (4,2 _ y 2 -wave SC and AF orders, and compared the re- 
sults with that in the dominant d 2 ,2_ a 2-wave phase with 
a higher harmonics term. In the coexisting phase, we 
found that the spin response is commensurate in a sub- 
stantial frequency range below a crossover frequency ui c 
for all dopings considered, and it switches to be incom- 
mensurate when the frequency is higher than lj c . The 
theoretical calculations are shown to be in good agree- 
ment with the experimental measurements. However, in 
the dominant d x 2_ y 2--wave phase with a higher harmon- 
ics term, the dispersion is just like that of the hole-doped 
one, namely exhibits a hourglass-like dispersion, which is 
not consistent with experiments. Thus, our result sug- 
gests that the inclusion of the two-band effect is impor- 
tant to consistently account for both the dispersion of 
the spin response and the non-monotonic gap behavior 
in the electron-doped cuprates. 
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FIG. 1: (Color online) Mean-field phase diagram for t — t' — 
t" — J model, where Ai is the superconducting order parame- 
ter without considering the AF order. The model parameters 
are taken as: t = -3.0J, t' = 1.02J, t" = -0.51J. 
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FIG. 5: (Color online) Frequency dependence of the intra- 
and inter-band contributions to the bare spin susceptibility 
Xo(q, w) [(al) and (a2) denote the imaginary part, (bl) and 
(b2) the real part] in the coexistence of AF and SC state at 
doping x=0.15. (al) and (bl) show the intra-band contribu- 
tions , and (a2) and (b2) the inter-band contribution. 



